1 Read system data and assign an initial approx. solution, or use flat-start 



2 Form Nodal Admittance Matrix, Initialize iteration count r=l 
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Compute Qp^'*'^by relation (6) after adjusting complex 
voltage for specified value by relation (10) 





no 

r 


Set SHp^PSHp+jQp^^*'^ 
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Compute V p from relation (5), and adjust complex 
voltage for specified value by relation (10) 



Compute V p^'**^ from relation (5) 
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Compute changes in the real and imaginary parts of the node-p 
voltage by using relations (7) and (8), and replace current value of 
DVRMX and DVIMX respectively in case any of them is larger 
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Calculate accelerated value of V p^^*'^ by relation 
(9) and update voltage for immediate use in the 
next node- voltage calculation by V p' = Vp^'"*^'^ 
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p=p+l 




Q STOP ^ 



Calculate line flows, and 
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and reactive power 
generation at PV-nodes 
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Fig* la: Prior Art: flow-chart of Gauss-Seidel Loadflow (GSL) Method 
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Read system data and assign an initial approx. solution, or use slack-start 



Form Nodal Admittance Matrix, calculate CosOp & SincDp using (19) and (20) & Initialize lPRP:=ITRQ=n=0 




YES 



CosOp=Cos (-48) & Sin(Dp=Sin (-48) 



Form admittance matrices [YG] & [YV], and factorize them & store as stated in calculation step-d 



Compute APj' & AO/, where i=l .... (m+V\ and i=l ....m 




8 Solve eqn. (13) for [AQ] and update voltage angles using, [G] = [G] + [AG] 



n 



Set voltage magnitudes of PV-nodes equal to the specified values, 
and hicrement ITRP=ITRP+l and r=(ITRP+ITR0)/2 



i Compute APj'& ACJi". where i=l,..,(nrH-k) and j=l,..,m 




Calculate line flows, and 
output indications of 
transformer tap position and 
reactive power generation at 
PV-nodes 



YES 



STOP ) 



j Compute FRO! using eon. (18) for only PO-nodes 



j_ 



Z 



Solve eqn. (15) for [AV] and update PQ-node magnitudes using [V] = [V] + [AV] 
While solving (IS), skip all the rows and columns corresponding to PV-nodes 



Adjust loadflow solution, if generator reactive power generation and transformer tap position limits are violated 



m Increment ITRP=ITRP+1 and r=(ITRP+ITRQ)/2 



Fig.lb:Prior Art: Flow-chart: Super Super Decoupled Loadflow (SSDL) Method 
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Fig. 2a: One-Une diagram of IEEE14-node network 




Fig. 2ci Non-redundant Level-1 sub-networks of fig. 2b are regrouped to reduce 
the number of processors required without increasing the number of 
nodes in any regrouped sub-network larger than the original largest 
sub-network of 6-nodes 
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Fig. 2b: Level-l sub-networks around circled nodes for the network 
of fig. 2 a 
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Read system data and assign an initial approx. solution, or use flat-start 



Form Nodal Admittance Matrix, Initialize iteration count r=l 



Initialize node count p=l , and maximum real and imaginary corrections DEMX=0.0 and DFMX=0.0 < 




yes 



yes 



Compute Qp^^^^^by relation (6) after adjusting complex 
voltage for specified value by relation (10) 




SetSHp = PSHp+jQ;'^**> 



SetSHp = PSHp+jQp"* 



SetSHp = PSHp +jQp'" 



Compute V p^'*'^from eqns. (27), (28), and (29) and adjust 
complex voltage for specified value by relation (10) 



Compute V p^'*" from eqns. (27), (28). and (29) ^ 
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Compute changes in the real and imaginary parts of the node-p 
voltage by using relations (7) and (8), and replace current value of 
DVRMX and DVIMX respectively in case any of them is larger 



99 
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Calculate accelerated value of V p^'*'^ by relation 
(9) and update voltage for immediate use in the 
next node-voltage calculation by Vp'= V p*'*** 



> n=r+l 



p=p+l 




Q STOP J 



Calculate line flows, and 

output indications of 
transformer tap position 

and reactive power 
generation at PV-nodes 
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Fig. 3a: invention: flow-chart of Gauss-Seidel-Patel Loadflow (GSPL) Method 
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A-PARALLEL-PROCESSOR 
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Read system data, assign an initial approx. 
solution, or use flat-start to initialize vector 
V(I), and initialize iteration count r=l 




Form Nodal Adniittance Matrix for a sub- 
network, and store it in the local memory 
of an assigned processor 




► 



Initialize maximum real and imaginary 
corrections DEMX=0.0 and DFMX=0.0, and 
VV(I)- CMPLEX (0.0, 0.0) 



read relevant elements of V(I) from 
commonly shared memory by all processors 




Perform two GSPL iterations for each sub- 
network. Steps-4 to 88 or 99 in Fig. 3a to 
calculate new estimate of voltage for each 
sub-network nodes 
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Calculate [VV(I)-V(I)], and Compute changes in the real 

and imaginary parts of the node-p voltage by using 
relations (7) and (8), and replace current value of DEMX 
and DFMX respectively in case any of them is larger. 



Calculate accelerated value of V V(I) by relation 
(9) and Perform operation V(1)=W(I) 



r=r+l 




29 



30 



Calculate line flows, and 

output indications of 
transformer tap position 

and reactive power 
generation at PV-nodes 



STOP ^ 



Fig. 3b J Invention: flow-chart of Parallel-Gauss-Seidel-Patel Loadflow 
(PGSPL) Method 
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Main/Server 
Processor, 
which is a 
delegating 
and 

coordinating 
processor 



< ► 
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< ► 
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Processor-n 
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Fig. 4: Invented Parallel computer Architecture/organization 
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Obtain on-line data 
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Initial Given/Established Specification of Controlled Parameters (Real 
and Reactive Power at PQ-nodes, Real Power and Voltage Magnitudes 
at PV-nodes, and Transformer turns ratio etc.) 



30 



Loadflow 
Calculation 



40 



Estimation 
of 

Power System State 



50 



NO 



60 




Evaluation, Good? 
(No over load & No 
under/over Voltage?) 




YES 



Correct/change One or 
More Controlled 
Parameters 



70 



Effect Actual Control 
To set finally obtained 
Controlled Parameters 



( END ) 



Fig. 5: Loadflow Calculation in Power Flow Control and/or 
Voltage Control in Electrical Power System 
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44 



Generator Loads 



Adjust Node 

MW 
Distribution 
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Assign Initial Approx. or Guess 
Loadflow solution, and Read data of 
Operating Limits of various Network 
Components, including Generator 
Excitation and Transformer Tap- 
changing Limits 



Machine 
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Perform Loadflow 
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& Solving any of 
the Load Flow 
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Load Flow Solution including 
Generators Reactive Power and 
Transformers Turns Ratio 
Indications 



Fig. 6: Load-Flow Calculation for Voltage Control in Electrical 

Power System 
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Fig* 7: An Exemplary 6-node Power System 

Nodes: 1, 2, 3 are PQ-nodes 

Nodes: 4 and 5 are PV-nodes 

Nodes: 6 is the slack/swing/reference node 

Transformers Tl and T2 are tap-changing 
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